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The kinetics of film growth of hydrates of methane, ethane, and methane-ethane mixtures were studied by 
exposing a single gas bubble to water. The morphologies, lateral growth rates, and thicknesses of the hydrate 
films were measured for various gas compositions and degrees of subcooling. A variety of hydrate film 
textures was revealed. The kinetics of two-dimensional film growth was inferred from the lateral growth rate 
and initial thickness of the hydrate film. A clear relationship between the morphology and film growth 
kinetics was observed. The shape of the hydrate crystals was found to favour heat or mass transfer and favour 
further growth of the hydrate film. The quantitative results on the kinetics of film growth showed that for a 
given degree of subcooling, the initial film thicknesses of the double hydrates were larger than that of pure 
methane or ethane hydrate, whereas the thickest hydrate film and the lowest lateral growth rate occurred 
when the methane mole fraction was approximately 0.6. 

Natural gas hydrates are a potential energy resource with huge reserves\ Gas hydrates can also find use in 
gas storage^"^, transportation^"^, and gas mixture separation^. However, gas hydrates may pose a plugging 
hazard during the production and transportation of crude oil and natural gas. Understanding the kinetic 
rate of gas hydrate formation is fundamental to both the utilisation of gas hydrates and the prevention of gas 
hydrate formation. In general, gas hydrate formation initiates at the water/guest- fluid interface. Once a hydrate 
crystal forms at the interface, it can grow along the lateral and normal directions to the interface, thereby forming 
a hydrate film with a certain thickness between the water and guest fluid^'^°. 

The morphology and growth rate of the hydrate film at the water/guest-fluid interface are significant to the 
overall study of hydrate formation kinetics. The morphologies of hydrates formed with a methane/ethane/ 
propane gas mixture in bulk liquid water and on the surface of a water droplet^^ were previously studied. In 
the study on hydrate crystal growth in bulk liquid water, it was found that the morphology of the hydrate crystals 
changed with the degree of subcooling and the methane mole fraction in the gas mixture^ \ In the study on hydrate 
crystal growth on the surface of a water droplet, it was found that the crystal size decreased with increasing 
methane concentration in the gas mixtures. In addition, the crystal size was smaller than that of simple methane 
hydrates^^. The macroscopic crystal growth morphology for hydrates formed from a 89.4% methane/ 10.6% 
ethane mixture was also studied with an emphasis on the investigation of induction time for hydrate formation 
on water droplets It was observed that as the methane mole fraction changed by no more than 10%, the 
influence of gas composition on the hydrate morphology could not be clearly distinguished. More importantly, 
these studies focused on the morphology of hydrate crystals, while the relationship between morphology and film 
growth kinetics was not addressed. 

The lateral growth of hydrate films at the gas/liquid or gas/solid interface has been studied experimentally and/ 
or theoretically by several groups^'^^"^^. However, most reported studies are for simple hydrates formed from a 
single guest. The initial thickness and lateral growth rate of the hydrate film are the key quantities to be measured 
to obtain the volumetric hydrate growth rate, which is often very difficult to measure. Selected measurements for 
the thickening rate of hydrate films in the normal direction to the interface after initial lateral growth exist^°"^^. 
The reported techniques for measuring film thickening rate include micrometry, microscopy^^"^^, laser inter- 
ferometry^^, and the magnetic resonance imaging (MRI)^^. Neither the laser interferometry nor MRI can deter- 
mine the absolute initial thickness of the hydrate film, although they can determine the thickening rate to some 
extent. In our previous work^^, we developed a method to measure the initial thickness of a hydrate film formed on 
the surface of a gas bubble exposed in water using a microscopic camera. A series of initial thickness data for 
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methane hydrate films were measured using this method^^, and it was 
shown that the initial thickness of the hydrate film was inversely 
proportional to the degree of subcooling^^. This finding is critical 
for modelling the lateral growth rate of the hydrate film based on 
the heat transfer control model proposed by Mori^^. 

In this study, the film growth kinetics of simple ethane hydrates 
and methane- ethane double hydrates were visually studied, includ- 
ing the dependence of the morphology, lateral growth rate, and 
thickness of the hydrate film on the gas composition and degree of 
subcooling. This study aims to establish a relationship between 
hydrate morphology and film growth kinetics. 

Results 

The reported results are based on changes in the gas composition and 
degree of subcooling (ATsub)- The degree of subcooling is defined as 
^^sub ~ T^"^ ~ T^""^, where T'^^'p is the experimental temperature and 
T^*^ is the hydrate formation equilibrium temperature at the experi- 
mental pressure (calculated with the Chen-Guo hydrate modeP^). 
Generally, differences in the morphology of hydrate films formed 
by different gases are clear at lower degrees of subcooling but dis- 
appear at higher degrees of subcooling. Figure 1 and Figures S3 to S6 
in the Supporting Material display diverse and texture-rich morphol- 
ogies of hydrate films formed from eight groups of different gases at 
three different degrees of subcooling. The images shown in the fig- 
ures are views of the hydrate films. All of the images were captured 
during lateral hydrate film growth or just after complete coverage of 
the gas bubble surface by the hydrate film. 

According to Figure 1, it is very interesting that the morphologies 
of hydrate films formed by gas mixtures with methane mole fractions 
between 0.190 and 0.620 are similar to that of pure ethane hydrate 
films but are quite different from that of pure methane hydrate films. 
This observation suggests that for these systems, nucleation and 
growth of the hydrate crystal is dominated by ethane. The texture 
of hydrate films formed at lower degree of subcooling (ATsub < 1 K) 
is clear, with individual crystals commonly appearing as long and 
stratified features (for more details, see videos 1 to 3 in the 
Supporting Material). At higher degrees of subcooling (2 K < 
A^sub < 3 K), the hydrate crystals are shorter and grain -like in shape 
(for more details, see videos 4 and 5 in the Supporting Material). As 
observed earlier, the size of the hydrate crystals decreases with 
increasing subcooling. This phenomenon is due to an increased rate 
of formation of new crystal seeds at higher degrees of subcooling. The 
emergence of new crystal seeds on the surface of a growing crystal 
hinders continued growth of a few large crystals at the expense of 
many small crystals^^. Increasing the degree of subcooling therefore 
shortens the growth period of each crystal. 

When the methane mole fraction of gas mixtures increases to 
0.834 and even 0.923, the crystal morphology becomes leaf-like or 
dendritic (for more details, see videos 6 and 7 in the Supporting 
Material), which is significantly different from that of hydrate films 
formed from pure ethane, pure methane, and mixtures with the 
methane mole fraction ranging from 0.190 to 0.620. This difference 
may be attributed to the formation of structure II (sll) me- 
thane-ethane double hydrates^^"^°. The structural transition of meth- 
ane-ethane double hydrates from structure I (si) to sll has been 
demonstrated by Raman spectroscopy^^ and in situ neutron diffrac- 
1[q^28,32 ^ methane composition range {ycud of 0.72 to 0.75. 
Conversely, the transition from sll to si occurs over a 7ch4 range 
of 0.992 to 0.994^^'^°'^^. The leaf-like morphology observed here is 
very similar to that of a sll hydrate film formed at a higher degree of 
subcooling (ATs^b ~ 3 K) on the surface of the water droplet from a 
methane- ethane-propane (90 : 7 : 3 molar ratio) gas mixture^^. The 
similarity of the results suggests that sll hydrates are being formed. 
Further determination of the structures of hydrate crystals would 
request Raman spectroscopy data in addition to the visual observa- 
tions in this work. 



It is very interesting to observe that the hydrate crystals formed at 
low degrees of subcooling (ATs^b < 1 K) always tend towards shapes 
favorable for heat or mass transfer and further fast hydrate film 
growth. The formation of stratified crystals with peaks and valleys 
certainly increases the total interfacial area between the hydrate and 
water, and maximising heat transfer. The dendritic shape of the 
hydrate crystals facilitates both heat transfer and mass transfer. 
Mass transfer resistance may become a critical problem for gas mix- 
tures with lower ethane concentrations. It has been shown that 
ethane can be easily removed from the gas phase and enriched in 
the hydrate phase when a gas mixture with a lower ethane concen- 
tration forms a hydrate^^'^^. The rapid formation and growth of 
hydrate crystals may reduce the local concentration of ethane. In a 
constant volume system, this reduction in local concentration would 
shift the thermodynamic equilibrium for hydrate formation, result- 
ing in decreased subcooling compared to the bulk system and a 
reduced hydrate growth rate. In addition, the branches of dendritic 
crystals increase the exposed areas to higher concentrations of 
ethane, resulting in continued crystal growth. Dendritic growth is 
widely observed in crystal growth^^'^^, and numerous explanations 
have been proposed to explain its occurrence^^. As such, the dendritic 
morphology of the hydrate crystals observed herein further adds to 
the intriguing field of dendritic growth theory. 

Another characteristic of hydrate films formed by gas mixtures 
with high methane mole fractions (0.834 and 0.923) is the hetero- 
geneous shapes observed. We find both small gravel-like crystals and 
larger dendritic or leaf-like crystals in the main body of the film. It is 
quite interesting to note the near perfect polyhedral single crystals 
among the dendritic crystals formed by gas mixtures with the highest 
methane mole fraction of 0.923 (Figure 2). The edges of these poly- 
hedral crystals are very clear and similar to single crystals of pure 
methane hydrate (white circled zone in Figure 2), suggesting they do 
not originate from the main dendritic crystals. (Instead, the polyhed- 
ral crystals in panels (e)-(h) in Figure 2 may originate from tiny 
crystals broken from the film crust). Hence, this phenomenon indi- 
cates the coexistence of si and sll hydrates. The formation of si 
hydrate is dominated by methane, while the formation of sll hydrate 
is dominated by both ethane and methane. The concomitant forma- 
tion of si and sll hydrates at the initial stage (the formation time is the 
same as the kinetic growth of hydrates at the gas bubble surface in 
Figure 2) has been observed by neutron diffraction and Raman spec- 
troscopy for methane-ethane gas mixtures with methane mole frac- 
tions of 0.95 and 0.93^^'^\ which is comparable to the gas composition 
in this study (methane mole fraction of 0.92). Our previous work on 
the Raman spectra of methane- ethane-tetrahydrofuran hydrates 
directly confirmed the coexistence of si and slP^. The coexistence 
of si and sll can also be inferred from the macroscopic partitioning of 
methane and ethane between the hydrate and vapour phases for 
systems with tetrahydrofuran in the liquid phase^^'^^. The appearance 
of methane hydrate crystals may be caused by the very high local 
concentration of methane near the main dendritic crystals due to the 
relatively faster consumption of ethane, as stated above. 

The lateral growth rate (vy) and initial thickness (^) of the hydrate 
film measured at 278.0 K are shown in Tables SI and S2 (see the 
Supporting Material) and are plotted in Figures 3 and 4. The lateral 
growth rate and initial thickness are calculated from the correlations 
proposed by Peng et al.^^ The lateral growth rate is correlated with 
Vf = \l/AT^, where i/^ is a lateral growth rate parameter and AT is the 
degree of subcooling. For the initial hydrate film thickness, the cor- 
relation is3 = k/AT, where k is the inversely proportional coefficient. 
As shown in Figures 3 and 4, there is very good agreement between 
the experimental data and correlations. 

Discussion 

The various hydrate film morphologies observed in this study indi- 
cate that the kinetic behaviour of the hydrate film depends on the gas 
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(a)l.OK 



(b) 1.2 K 




(c) 1.3 K 



(d) 1.5 K 




(g)1.5K 



(h) 2.5 K 



Figure 2 | The single crystals in the growth of hydrate films formed with 
methane-ethane hydrate forming gases with methane mole fraction of 

0.923. Images (a)-(d) were recorded in the film lateral growth process and 
images (e)-(h) in the film regrowth process in the crust gap that was 
originated by supplying more gas to break the film crust after the gas 
bubble was fully covered by the hydrate film. The scale bars correspond to 
1 mm. 

composition and degree of subcooling. For example, there is little 
difference between the lateral growth rate data for methane hydrates 
in this work and the results obtained in our previous work using the 
same experimental method^^. These growth rates are only slightly 
higher than those obtained at a planar water/methane interface^^'^". 
However, these values are much lower than the lateral growth rates of 
methane hydrate films growing at the interface between methane and 
a sessile water droplet^\ The differences in lateral growth rate for 
different systems (gas bubble or water droplet) should be caused by 
different volumetric growth rates, heat production rates, and heat 
transfer environments. For hydrate film growth at a given degree of 
subcooling (AT) in this work, the lateral growth rate and initial 
hydrate film thickness can be readily estimated using ij/ and /c, 
respectively, for a known composition of the hydrate forming gas. 
The values ij/ and k for each hydrate-forming gas mixture are 
regressed from the experimental data in Figures 3 and 4 and plotted 
in Figure 5 (actual values are given in Tables S3 and S4 in the 
Supporting Material). As shown in Figure 5, the dependence of i/^ 




Figure 3 | Variation of lateral growth rate of hydrate film with subcooling 
for different hydrate forming gases. Solid Knes denote the results 
calculated from the correlation Vf = ij/AT^. 

and k on the gas composition is contrasting: the lateral growth rate of 
simple methane or ethane hydrate is higher than that of methane- 
ethane double hydrates, and the initial film thickness of simple meth- 
ane or ethane hydrate is smaller than that of methane- ethane double 
hydrates for similar degree of subcooling. The maximum 3 or min- 
imum Vf occurs at a methane concentration of approximately 0.62. 
These results lead to the conclusion that a lower lateral growth rate 
leads to a thicker hydrate film, as more time is available for vertical 
development of the film front. Another way to interpret this result is 
that a thicker hydrate film is unfavourable for heat transfer and 
therefore lowers the lateral growth rate. 

The product of 3 and provides the volumetric growth rate (r) of 
the hydrate at the gas -water-hydrate contact line. 



r = il/kATi = il/'Af2 



(1) 



where ij/' is the variation of the formation rate constant. The depend- 
ence of i/^' on the composition of the hydrate forming gas is shown in 
Figure 6 along with the volumetric hydrate growth rate. From the 
data, one can conclude that the volumetric growth rate of double 
hydrates is also lower than that of simple methane or ethane hydrate 
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Figure 4 | Initial hydrate film thickness measured for different degrees of 
subcooling and hydrate forming gases. Solid Hues denote the results 
calculated from the correlation 3 = k/AT. 
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— □ — parameters ^/ 
— ■— parameters A 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Methane content in the hydrate forming gas 

Figure 5 | Values determined for the parameters and ic for the different 
compositions of the hydrate forming gas. 



at a given degree of subcooling. However, the influence of gas com- 
position on the volumetric growth rate is not as pronounced as the 
influence of gas composition on the lateral growth rate. As discussed 
previously, the gas composition also affects the morphology and 
consequently the porosity of the hydrate film. Therefore, r is actually 
the apparent growth rate of the porous hydrate film and not the true 
volumetric growth rate of the hydrate lattice. 

From the analysis of the growth rate and thickness of the hydrate 
film, it can be observed that the lateral hydrate film growth rate is 
determined by both the seeding rate of crystal seeds and the growth 
rate of the single crystal, which were determined from these mor- 
phology and kinetic studies. Generally, the seeding rate of crystal 
seeds is very low at low degrees of subcoolings so the lateral growth 
rate is controlled by the growth rate of a single crystal. In this case, the 
low lateral growth rate leads to a thicker hydrate film. Conversely, the 
seeding rate of crystal seeds is high at high degrees of subcooling and 
the lateral growth rate is controlled by both the seeding rate of crystal 
seeds and the growth rate of single crystal. This results in a higher 
lateral growth rate and a thinner hydrate film. In addition to the 
degree of subcooling, the gas composition also has a strong effect 
on the seeding rate of the crystal seeds. At a given degree of subcool- 
ing, three observations can be made: i) methane hydrates show the 
highest seeding rate; ii) the seeding rate decreases with increasing 
methane mole fraction {ycH4 < 0.62), suggesting that hydrate nuc- 
leation is controlled by ethane in these cases; and iii) the seeding rate 
increases with increasing methane mole fraction (ycH4 > 0.62), sug- 
gesting that hydrate nucleation is controlled by methane in these 
cases. 

As stated previously, the gas composition dependence of the film 
morphology also indicates that the nucleation of hydrate crystal is 
dominated by ethane in the range 7ch4 < 0.62, while it is dominated 
by methane or methane and ethane together in the range ycH4 > 
0.62. In the range 7ch4 < 0.62, the driving force for the formation of 
si pure ethane hydrate crystal decreases with decreasing ethane con- 
centration. Additionally, hydrate formation results in a local concen- 
tration of ethane near the film front lower than that of the bulk gas 
phase. This effect becomes more obvious with decreasing ethane 
concentration, as stated previously. These two factors may lead to 
a decreased seeding rate of the crystal seeds with decreasing ethane 
concentration in this concentration range. However, in the range 
7cH4 > 0.62, sll hydrate was assumed to be formed. In this case, as 
there are more small cavities in sll than in si, methane likely has a 
greater contribution to hydrate nucleation, and therefore, the seeding 
rate increases with increasing methane concentration. These factors 
may explain why the lowest lateral growth rate occurs around )/ch4 = 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Methane content in the hydrate forming gas 

Figure 6 | Variation of formation rate constant (^0 with composition of 
the hydrate forming gas. 

0.62. The lowest lateral growth rate certainly corresponds to the 
largest thickness of the hydrate film. 

In summary, the hydrate film growth processes of pure methane 
hydrate, pure ethane hydrate, and six methane- ethane double 
hydrates at the gas bubble/ water interface were studied using micro- 
scopic visual methods. The morphologies, lateral growth rates, and 
thicknesses of the hydrate films were measured as a function of the 
gas composition and degree of subcooling. Based on these measure- 
ments, insight into the morphology and growth kinetic behaviour of 
methane -ethane hydrate films can be summarised as follows: 

• (i) For )/cH4 = 0.0 to 0.62, the morphologies of the hydrate films 
formed from gas mixtures are similar to that of simple ethane and 
quite different from that of pure methane. This result suggests 
that the nucleation and growth of the hydrate crystals is domi- 
nated by ethane. The observed crystals change from stratified to 
wheat ear-like with increasing subcooling. 

• (ii) For 7cH4 — 0.834 to 0.923, the crystal morphology becomes 
leaf-like or dendritic at lower degrees of subcooling, indicating the 
formation of sll hydrate. Additionally, polyhedral crystals like in 
pure methane hydrate can be observed among the main dendritic 
hydrate crystals, indicating the coexistence of si and sll hydrates. 

• (iii) The hydrate crystals formed at low degrees of subcooling 
always tend towards shapes that minimise heat or mass transfer 
and further favour hydrate film growth. 

• (iv) The lateral hydrate film growth rate is controlled by the 
seeding rate of crystal seeds and the growth rate of single crystals. 
Both the degree of subcooling and gas composition strongly affect 
the seeding rate and therefore strongly affect the lateral hydrate 
film growth rate. At lower degrees of subcooling, the seeding rate 
of crystal seeds is very low and film growth is controlled by the 
growth of single crystals. At higher degrees of subcooling, the 
seeding rate of crystal seeds is high, and the lateral growth rate 
is controlled by both the seeding rate of crystal seeds and the 
growth rate of a single crystal. At a specified degree of subcooling, 
methane hydrates show the highest seeding rate. The seeding rate 
decreases with the decreasing ethane mole fraction in the range of 
7cH4 from 0.0 to 0.62, which further indicates the nucleation of 
hydrate is controlled by ethane in these cases. For higher 7ch4 
(> 0.62), the seeding rate increases with decreasing ethane mole 
fraction, suggesting that hydrate nucleation is controlled by 
methane in these cases. 

• (v) Lateral film growth rate and initial film thickness can be 
correlated by a kinetic model based on the assumption that heat 
transfer is the controlling step. The initial hydrate film thickness 
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and lateral hydrate film growth rate are found to vary inversely 
with the gas composition. 



Methods 

The experimental apparatus and procedure are the same as those used in our previous 
study^^'^^; details are provided in the Supporting Material. Double-distillated water 
was used in all measurements. Methane and ethane, both with a purity of 99.9% (mole 
basis), were provided by Beijing Beifen Gas Industry Corp. Six methane + ethane gas 
mixtures, with methane mole fractions of 0.190, 0.350, 0.519, 0.620, 0.834, and 0.923, 
were synthesized in our laboratory (mixed based on partial pressures) and their 
compositions were analysed by gas chromatography (HP6890). For all experiments 
performed, the temperature of the system was set to 278.0 K. The pressure for each 
experiment was varied to achieve the desired degree of subcooling. 
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